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Abstract 

We review the principle and methodology of leakage radiation microscopy (LRM) applied to 
surface plasmon polaritons (SPPs). Therefore we first analyse in detail the electromagnetic theory 
of leaky SPP waves. We show that LRM is a versatile optical far-field method allowing direct 
quantitative imaging and analysis of SPP propagation on thin metal films. We illustrate the LRM 
potentiality by analyzing the propagation of SPP waves interacting with several two dimensional 
plasmonic devices realized and studied in the recent years. 

PACS numbers: 
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I. INTRODUCTION 



In recent years intensive investigations of surface plasmon polaritons (SPPs) have been 
made in the promising context of nanophotonics. This research is actually motivated by the 
current trends for optical device miniaturisation and by the possibilities of merging aspects 
of nanophotonics with those of electronics. SPPs are electromagnetic waves bounded 
to dielectric-metal interface. As surface waves, SPPs are exponentially damped in the 
directions perpendicular to the interface Furthermore, SPPs could be used to transfer 
optical information in a two dimensional (2D) environment. This appealing property 
can be used for optical addressing of different 2D optical systems and nanostructures 
located at a dielectric/metal interface. Actually several 2D SPP devices including passive 
nanostructures including mirrors or beam splitter and active elements like molecules or 
quantum dots are currently under development and investigation. Developments such as 
these raise the prospect of a new branch of photonics using SPPs, for which the term 
"plasmonics" emerged {2 4]. 

However, for experimental investigations of optical devices an important characteristic 
of SPP modes is that their spatial extent is governed and defined by the geometry of the 
nanoelements rather that by the optical wavelength [7| . This consequently opens possibilities 
for breaking the diffraction limit but requires instruments of observation adapted essentially 
to the subwalength regime and being capable of imaging the propagation of SPPs in their 2D 
environment. Usually the analysis of the subwalength regime implies necessarily near field 
optical (NFO) methods {5, 6] able to collect the evanescent (i. e., non radiative) components 
of the electromagnetic fields associated with SPPs. However, when the metal film on which 
the 2D optical elements are built is thin enough (i. e ., with a thickness below 80-100 
nm) and when the subtratum optical constant (usually glass) is higher than the one of 
the superstratum medium an other possibility for analyzing SPP propagation occurs. This 
possibility is based on the detection of coherent leaking of SPPs through the substratum. 
Such a far-field optical method is called leakage radiation microscopy (LRM) fl-lio] and 
allows indeed a direct quantitative imaging and analysis of SPP propagation on thin metal 
films. 

The aim of this article is to present a short overview of recent progress in the field of SPP 
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imaging using LRM. In a first part of this work we will describe the theoretical principles 
underlying LRM. In the second part we will discuss modern leakage radiation methods and 
illustrate the LRM potentialities by analyzing few experiments with SPP waves interacting 
with 2D plasmonic devices. 



II. LEAKAGE RADIATION AND SURFACE PLASMON POLARITONS 

In order to describe the theoretical mechanisms explaining leakage radiation it will be 
sufficient for the present purpose to limit our analysis to the case of a metal film of complex 
permittivity €i(oj) = e[ + ie" (u = 2ttc/X is the pulsation) sandwiched between two dielectric 
media of permittivity e (substrate) and e 2 < eo (superstratum). This system is theoretically 
simple and to a good extent experimentally accessible [l[ In the limiting case where the 
film thickness D is much bigger than the SPP penetration length in the metal (i. e., D > 70 
nm for gold or silver in the visible domain) one can treat the problem as two uncoupled 
single interfaces. We will consider as an example the interface 0/1 (the media and 1 are 
located in the domain z > and z < respectively ). Such an interface will be identified in 
the following with the plane z = const, in cartesian coordinates. An elementary harmonic 
SPP wave is actually a TM electromagnetic mode characterized by its pulsation u and its 
magnetic field H = [0, H y , 0] where the y component can be written 

H = ae ikxX e ikz0Z e~ iu)t in the medium 
H x = e *.* e <*.i* e -krt i n the medium 1, (1) 

and where k x = k' x + ik" is the (complex valued) wavevector of the SPP propagating in 
the x direction along the interface. k Z j = kj = ±A/[(w/c) 2 e i — are the wave vectors in 
the medium j =[0 (dielectric), 1 (metal)] along the direction z normal to the interface. By 
applying boundary conditions to Maxwell's equations one deduces additionally a = 1 and 

/Ci kn , . 

- - - = 0, 2 

ei e 

which implies 



k x = ±(u/c)J^- (3) 
V e o + e i 



k 3 =±{u/c)\—t- (4) 

fcQ T tl 



for a SPP wave propagating along the x direction. The choice of the sign convention con- 
necting the z and x components of the wave vector is a priori arbitrary and must be done 
only on a physical ground. Indeed, due to ohmic losses in the metal we expect an expo- 




FIG. 1: Structure of the SPP magnetic field Realty] across an interface air /gold (thick line) and 
glass/gold (dashed line). The optical wavelength considered is A = 800 nm. The permittivity of 
glass is taken to be e g i ass = 2.25. The inset shows the conventions for the axes x and z. The arrows 
indicate the direction of the real part of the wave vector normally and parallely to the interface. 

nentially decaying SPP wave propagating along the interface. This condition implies the 
relation k' x ■ k" x > [ll]. This inequality is actually always fulfilled since from Eq. 3 one 
deduces 

which is indeed positive because e'[ > 0. By writing k Z j = k'- + ik" one additionally obtains 
the relation 

- k' ■ kl = (u/c) 2 e';/2 - k[ ■ k'l = k' x ■ K > 0. (6) 

This relation fixes the sign conventions since the wave must also decay exponentially when 
going away from the interface in both media. More precisely one gets 

K ■ K < o, (7) 

K ■ K => 0. (8) 

The product k[ ■ k'[ is positive if e[ > — |ei| 2 /(2e ), a fact which is indeed true for silver and 
gold interfaces with air or glass in most of the visible optical domain. However small negative 



values of Eq. 8 occur for silver close to the interband region around A ~ 350nm. Additionally 
a higher value of eo will also change the sign in Eq. 8. Fig. 1 shows the behavior of the SPP 
magnetic field close to an interface gold/air and gold/glass at the optical wavelength A = 800 
nm. At such a wavelength the conditions given by Eqs. 5-8 impose the solutions 

k = ±(oVc)y/^_ k . = _ (w/c)v C]L (9) 

The real parts of the k z components of the SPP wave vector are for both media oriented in 
the same direction corresponding to a wave propagating from the air side to the metal side 
(see inset in Fig. 1). Furthermore the waves are exponentially damped when going away 
from the interface in agreement with Eq. 7, 8 (see Fig. 1). Most important for us is that the 
Poynting vector [12|] S = Real[E x H*]/2 is defined in the medium j by 

S,- = -cBjeal[ kxk+ ! ej± ]e- 2k ^- 2k " z . (10) 

On the dielectric side the energy flow is as expected oriented in the direction of Real[k]. 
However it can be shown on the metal side and for wavelengths not too close from the 
spectral region associated with the interband transition of gold or silver that the energy 
flow in the x direction is oriented oppositely to the wave vector RealffcJ since Real^/ei] = 
(k' x €[ + k' x e'{) / \ei\ 2 is dominated by k' x ■ e[ and since e[ < 0. However the total energy flow 
in the x direction S x = / + °° Sx dz + Sl metal W is oriented along k' x . Additionally in 
the z direction Real[&i/ei] = (A^e^ + k"e'[)/\ei\ 2 is parallel to k[ since k"e" dominate. This 
implies in particular that the energy associated with the SPP wave is absorbed by the metal 
during its propagation through the interface from the air side to the glass side. It should 
be observed that already in the case of the ideal plasma model neglecting losses the wave 
vector k x = k' x is antiparallel to S x in the metal since there is no imaginary part and since 



e x = 1 - u 2 /u 2 < (see also [13|, 



We show in Fig. 2 the curves associated with the dispersion relations of SPPs propagating 
along a gold/air and gold/glass interface respectively. Fig. 2 A represents the dependencies 
u versus k' x corresponding to Eq. 3. Fig. 2B shows the dependencies u versus L S pp where 
Lspp = 1/ (2A£) is the propagation length of the SPP waves (for the metal optical constant 



we used the experimental values given in 15|). The typical back bending around A 



520nm corresponds to the resonance associated with the bound electrons and the interband 



transition (for a good discussion see |17j). Far away from the interband the real part of 



the dispersion is close to the asymptotic light lines: we speak about Zenneck surface modes 
by opposition to Fano and evanescent modes existing close to the interband 13|, [l6j. An 




FIG. 2: Dispersion relations of SPP waves propagating along a air/gold (thin dashed curve) and 
glass/gold (thick dashed curve) interface. (A) Real part of the dispersion relations showing the 
evolution of k' x with uj/c = 2ir/\. The light cones corresponding to an optical wave propagating 
in air and glass (i. e., in the bulk medium) are represented by a thin and thick continuous line 
respectively. The back bending at A = 520 nm correspond to the interband resonance. (B) 
Imaginary part of the dispersion relations showing the evolution of L$pp with oo/c = 2tt/X. The 
damping is higher for the glass/gold interface. 

important feature occurs close to this interband transition since the slope in Fig. 2A is 
diverging. This means that the group velocity defined by v g = du/dk' x is infinite. One can 
qualitatively deduce that there is however no paradox with causality by remarking that the 
propagation length is approaching zero at such wavelengths. Actually the SPP life time 
rspp defined by v g ■ T$pp ~ Lspp is tending as well to zero. This is a clear signature of the 
absence of significant SPP propagation in the interband spectral region. This fact actually 
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precludes faster than light information transfer by SPP waves and is in agreement with 
relativistic causality. 

A second important feature concerning Fig. 2 A is that the air/gold dispersion curve is 
located inside of the light cone for glass defined by the equation K 2 = (w/c) 2 e g i ass where K 
is a real light wave vector and e g i ass = 2.25. Writing K = [K x , 0, K z ] the wavevector of a TM 
(i. e., p polarized) plane wave propagating away from the interfaces into the glass side one 
see that SPPs propagating at the air/metal interface can radiate into the glass substrate if 
the condition 

K x a±Real[(u/c)J-^-} (11) 
V e o + e i 

is approximately fulfilled. Here we neglected the role of the imaginary part in Eq. 3. Similarly 
one can deduce that none of the SPPs propagating at the two interfaces can radiates into 
the air side. 

In order to have a more complete analysis one must actually consider the problem with two 
coupled interface 0/1 (glass/metal) and 1/2 (metal/air) supporting SPP waves and separated 
by a small distance D. The two interfaces are coupled by evanescent SPP waves tunnelling 
through the metal slab. Such a mathematical problem can only be treated numerically by 
resolving an implicit equation. As for the single interface this equation can be defined directly 
from Maxwell Equations 111]- However it is much easier and convenient for the following 
to remark with Raether |l| that Eq. 2 and consequently Eq. 3 are obtained by finding the 
zeros of the numerator in the Fresnel reflectivity coefficient for a TM wave coming from the 
dielectric side: 

_ (k /e - fci/ei) , . 

^-(^Ai + Weo)' { } 

with kj = ±a/ ((u/c) 2 ej — A; 2 ,). Actually Raether |l| reasoned with the denominator of the 

Fresnel coefficient due to different conventions for the signs of the wave vectors kj. However 

it is remarkable that the result is the same at the end of the calculations. Identically one 



can thus define the Fresnel coefficient for a TM wave reflected by the slab 0-1-2 |l|, 1 121. Il7||: 



identically 

1-2 fl y 



°< 1 ' 2 " l + R p 0jl R p 1>2 e^ D ' [ ' 

and find the zeros of the numerator, i. e., one can solve the implicit equation 

R p 01 + R p 2 e 2iklD = 0, (14) 
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in order to define the SPP dispersion relations. From this equation it follows that 

(k /e + k 1 /e 1 )(k 2 /e 2 -k 1 /e 1 )e ikiD 
+(fc 1 /e 1 - k /e )(k 2 /e 2 + h/e^e' 1 ^ = 0. (15) 

As for the single interface one has an important relation between the real and imaginary 
parts of the SPP wave vectors in the different medium: 

- k' ■ K = ("/cfe'l/2 - k[ ■ k'l = -k' 2 ■ k" = k' x ■ k" x . (16) 

Since we are interested only into the solutions which are decaying along the interface we 
(in agreement with our previous treatment of the single interface) suppose the condition 
k' x ■ k" x > satisfied. A second important point is that due to the arbitrariness in the sign 
of kj there are in fact apriori 8 possibilities for writing Eq. 15. However, Eq. 15 is invariant 
under the transformation k\ — > —k\. This means that the number of apriori possibilities 
or the sign of kj is reduced from 8 to 4. This multiplicity was studied by Burke et al. 



ll| however since for the present purpose we are looking for SPP waves leaking from the 
air/metal interface into the glass substrate we will consider only the possibility 

k = -y/((u/c)*e -k*) 
h = +V(("/c) 2 e 1 - hi) 

k 2 = -^/{{uj/cYe 2 -kl). (17) 

The sign of k\ is however arbitrary as explained above and we choose it here positive by 
definition. In order to define a SPP wave leaking into the glass substrate one has thus 
;o substitute Eq. 17 into Eq.15 and find numerically (i. e., by a minimization procedure 
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19| ) the zeros of the implicit equation with variable k' x and k". This has to be done 
only in the quarter of the complex plane corresponding to k' x > 0, k" x > 0. The quarter 
k' x < 0, k" x < must be equivalent due to symmetry and corresponds actually to decaying 
SPP waves propagating in the negative x direction. The two other quarters of the complex 
plane correspond to growing SPP waves along the interface and will be rejected on a physical 



ground (compare 

Figs. 3A and 3B show numerical calculations of dispersion relations corresponding to a 
SPP wave leaking through a gold film from the air side to the glass side. The thickness 
is taken to be D=70, 50, 20, and 10 nm respectively. For the value D > 70 nm the 
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FIG. 3: Dispersion relations of SPP waves propagating along an air /gold interface and leaking 
into glass. The dispersion is calculated for four film thickness D : 70 nm (dashed line curve)es, 50 
nm (continuous line with triangular markers), 20 nm (continuous line with square markers), and 
10 nm (continuous line with circular markers). (A) Real part of the dispersion relations showing 
the evolution of k' x with uj/c = 2tt/X. The light cone corresponding to optical wave propagating 
in air (i. e., in the bulk medium) is represented by a continuous line. (B) Imaginary part of the 
dispersion relations showing the evolution of L$pp with uj/c = 2ir/\. 

dispersion relation is identical to the dispersion for the single air/gold interface for semi 
infinite media. However, for smaller thickness the coupling between the interface increases 
and the propagation length decreases as shown on Fig. 3 B. The magnetic field associated 
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FIG. 4: Structure of the magnetic field Real[-ff y ] associated with a SPP mode leaking through a 
gold film of thickness 50 nm (medium 1) from the air side (medium 2) to the glass side (medium 
0). The evolution is represented along the normal z to the interfaces. The optical wavelength 
considered is A = 800 nm. The two horizontal lines show the interfaces separated by 50 nm. The 
inset shows the conventions for the axis x and z. The arrows indicate the direction of the real part 
of the wave vector normally and parallely to the interface. 

with SPP electromagnetic modes in this layered system is given by 

JJ _ e ik x x e ik z g—ivt 

H x = e ikxX [asm (k lZ ) + /3cos (hz^e^ 
The coefficients a,/3,7 are obtained by considering the boundary conditions and one finds 

• foci a 1 

a = i- , p — 1 

fcie 

k e 

7 = z-^sin(A;iD)+cos(A; 1 D). (19) 

As an illustration we show in Fig. 4 the evolution of the real part of the magnetic field across 
a 50 nm thick gold slab sandwiched between the glass substratum and the air superstratum 
for an optical wavelength A = 800nm. As visible the SPP field is located in the vicinity of 
the air /gold interface and is evanescent on the air side. This is clearly reminiscent of our 
former analysis of SPPs propagating at the single air/metal. In addition however the wave 
is leaking radiatively (i. e., propagatively) into the glass substrate. However from Eq. 18 
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and the condition k' ■ k'l < it is clear that the leaking wave is exponentially growing in 
the -z direction when going away from the gold slab. This is already the result we obtained 
when we considered the limit of the thick slab. An exponentially growing wave looks non 
physical and is in particular associated with infinite radiated energy in the far field. There 
is now the question of how to connect a growing wave with the basic reasoning giving rise to 
Eq. 11 and the idea of phase matching between the (real part) of the SPP wave vector with 
a propagative plane wave vector in the glass substrate. However such paradoxes disappear 
for two reasons: First, an infinite energy occurs only because we considered an infinite 
interface or equivalently because we did not consider how the SPP is locally launched on 
;he metal film. When such conditions are taken into account this paradox must disappear 



111 ] . Second, the SPP wave defined by Eq. 18 is actually a wave packet when looked at 
through the Fourier basis of propagative TM plane waves. Since in the far field (i. e., in the 
glass substrate) one actually detect such plane waves one must do a Fourier transform in 



order to generalize Eq. 11 



Instead of Eq. 11 one obtains consequently a statistical 



distribution of (real) wavevectors K x given by 



const. . . 

'<*■> = (20) 

where 2k" defines the full width at half maximum (FWHM) of this Lorentz distribution of 
radiated plane waves. By noting as usual 9 the angle between the wave vectors K of the 
plane waves refracted into the glass substrate and the normal z to the interfaces one has by 
definition K x = 27r v /io/A sin# and the angular distribution of radiated power is in the far 
field given by: 

T (0] = (OA \ 

K ] (27i^/\sm9-k' x y + (k^y 1 ; 

III. LEAKAGE RADIATION MICROSCOPY 



Historically the first observations of leakage radiation by SPP propagating on a thin metal 
film were reported by analyzing scattering by rough metal surfaces of light into SPPs [l|, l^cl . 
The possibility of using rough surface to excite SPPs was extensively studied in the past [1] 
and is based on the fact that the scattering by small defects on a flat film can represent a 
source of evanescent momentum sufficient for the light waves to match the SPP dispersion 
relation. Equivalently the amount of momentum needed can be carried by grating coupling 
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SPP waves are subsequently emitted back into the glass substrate as leakage radiation 
(see Fig. 5). This light collected on a photographic plate forms a ring-like distribution 





FIG. 5: Observation of leakage radiation (LR) through a rough thin silver film (wavelength A = 
550nm). (A) SPPs are launched on the air /silver interface by using scattering by film rugosities 
to produce the amount of momentum necessary to match the SPP dispersion relation. SPPs are 
leaking through the glass substrate and detected using a photographic plate in the far-field. (B) 
Photography showing the le akag e radiation cone. The photography is extracted from the work by 
H. J Simon and J. K. Guha [2o| (© Opt. Comm., Elsevier, 1976). 



in agreement with Eqs. 11, 20, 21. The FWHM of the S 



'P wavevector distribution is in 



direct correspondence with the radial width of the ring 20]. Further progress was obtained 
recently with the development of near field scanning optical microscopy (NSOM) which 
allows the local optical excitation of evanescent waves in the vicinity of a metal surface. 
Such evanescent waves can carry a sufficient amount of momentum to couple to SPP waves. 



Direct observations have indeed confirmed this principle |8|,[21|, 22]. As an example we show 
in Fig. 6 an experiment in which the NSOM tip launches SPPs on an aluminum film which 
after interaction with a hole excites optically some quantum dots (QDs) located below. The 
collected signal shows a specific QD luminescence spectrum j^. By scanning the sample 
around the NSOM tip one can realize SPP mapping since the hole acts a probe structure 
for the field emitted by the tip. Quantitative analysis of the total luninescence of the QDs 
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associated with a given hole show clearly that the QD excitation is mediated by SPPs 
propagating on the aluminum film. Fig. 6B shows the radial dependence of the collected 




FIG. 6: (A) Sketch of the experiment using a NSOM tip to excite SPPs on an aluminum film 
with nanoholes. The film containing 240 nm diameter holes covers a CdTe/ZnTe system of QDs 
whose spectral luminescence induced by SPP excitation is characteristic of a given hole. The laser 
excitation wavelength is 514.5 nm. (B) Radial mapping of the SPP intensity. The SPP intensity 
is proportional to the QDs luminescence of a given hole. The curve represents the evolution of the 
total luminescence collected (using a detector in the optical far field) as a function of the distance 
R separating the NSOM (source) from the hole (probe). The experiment is made at 4.2 K [22| (© 
Europhys. Lett., EDP, 2003). 

intensity. These results agree well with a 2D SPP dipole model supposing an effective dipole 



ocated at the tip apex [8|, |22j (see also the discussions concerning t 
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26] theory of diffraction by a small aperture in a metal film in 



re Bethe-Bouwkamp 
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31]). Following this 



model the SPP wave can be theoretically modelled by a scalar wave \&(p, 6) given by 

— , „ e ifcsppp cos0 
(p, 0) = const. — — (22) 

VP 

where /cspp = k' x + ik" is given by Eq. 3, p, 4> are polar coordinates on the metal film and the 
origin of the coordinate is taken at the dipole position. is the angle between the dipole 
associated with the NSOM tip (parallel to the polarization of the laser beam injected in the 

13 



NSOM tip) and the the radial vector p = [x,y]. This simplified model can be theoretically 
justified by using the Green Dyadic Formalism {s^J and has been applied by many authors 



22 



33 



36[ to SPP waves propagating in various environments. 



n 



successfully 

Several authors applied NSOM methods coupled to LRM jg, 9]. In particular in [8| Hecht et 
al. realized an optical setup using an immersion oil objective to collect the leakage radiation 
emitted by the NSOM tip on gold or silver films (see fig. 7). The system shown in Fig. 7 B is a 



(A) 




NSOM 



10 pm 



FIG. 7: (A) Principle of the experiment to generate SPP leakage radiation (LR) using a NSOM 
tip. LR is collected using an immersion oil objective (not shown here). (B) Map of SPP intensity 
using a charge coupled device (CCD) camera to collect the LR emitted through the glass substrate. 
SPPs are launched on a 60 nm thick silver film at the optical laser wavelength A = 633nm. The 
image in B is taken from B. Hecht et al. [8] (© Phys. Rev. Lett, American Physical Society, 1996). 



60 nm thick silver film optically excited by a NSOM tip at the laser wavelength A = Q33nm. 
It can be shown by analyzing Fig. 7 B that the radiation pattern is well described by a 
2D dipole model in agreement with Eq. 22. In particular the SPP propagation length was 
measured and is in fair agreement with our analysis in Sect. 1. Additionally it was shown in 
js] that one can also analyze the Fourier distribution of SPP momentum (given by Eqs. 20, 
21) by defocusing the objective lens. As expected SPP rings similar to the one of Fig. 5 
were observed. 
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In the same context we developed in recent years a systematic approach using far-field 




FIG. 8: (A) Principle of the experiment to generate SPP leakage radiation (LR) using a focussed 
laser beam (objective 0\, lOx, numerical aperture NA = 0.3). LR is collected using an immersion 
oil objective (objective O2, 63x, NA = 1.3) and refocussed on the CCD camera by using an auxiliary 
lens. (B) Map of SPP intensity using a CCD camera to collect the LR emitted through the glass 
substrate. SPPs are launched on a 50 nm gold film at the optical laser wavelength A = 800nm 
from a gold ridge (50 nm height, 150 nm width) represented by a white line. The size of the laser 
spot with a diameter ~ 8 /xm is represented by the white dashed circle. 

microscopy to analyze quantitatively the interaction between SPPs and plasmonic devices 
by using LRM. The nanodevices studied were all fabricated by electron beam lithography 
(EBL) allowing the precise and reproducible tailoring of metal and dielectric surfaces on a 
lateral size dimension down to 20 nm [3^]. As an example we show in Fig. 8 B a LRM image 
obtained by using a gold ridge (50 nm height, 150 nm width) lithographed on a 50 nm thick 
gold film to launch two well collimated and counter propagating SPP beams. These beams 
are launched by focussing a laser beam with a microscope objective (lOx, numerical aperture 
NA = 0.3) on the gold ridge. Scattering by the nanostructure gives rise to evanescent waves 
supplying the right amount of momentum necessary for generating a SPP wave. The optical 
LRM setup is sketched in Fig. 8 A. Leakage radiation emitted through the glass substrate 
is collected by an immersion oil objective (63x, numerical aperture, NA = 1.3). Light is 
subsequently refocussed on a charge coupled device (CCD) camera. The direct mapping of 
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the SPP intensity with this method provides a one-to-one correspondence between the 2D 
SPP intensity and the image recorded on the CCD camera. It should be observed that an 
incident laser beam with diameter Dx is in the focal plane (object plane) of the objective 
0\ focussed into a disc of diameter w (i. e., w =beam waist) such that 

2A Dx 

tana = — = — - (23) 

where / is the focal length of the objective and a the divergence angle of the laser beam 
focussed on the sample. The direct application of this version of the Heisenberg relation 38 1 
implies that the divergence angle otspp of the SPP beam launched on the metal film must 
equal the divergence angle a of the impinging laser beam. This result is in good agreement 
with the experimental case shown in Fig. 8 with Dx = 2 mm, w = 8.3 fim and a = 3.5°. 
Changing the objective focal length is a straightforward means to obtain different divergence 
angles aspp (see for example 
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39|,|40|). 



As a further improvement it is possible to modify the previous optical setup in order to 




FIG. 9: Experimental scheme for LRM. SPPs are excited by laser light focused by the microscope 
objective 01 (50x, numerical aperture, NA = 0.7) onto a structured gold thin film on a glass 
substrate. LR emitted into the glass substrate from the gold/glass interface £ is collected by the 
immersion microscopy objective 02 (63x, NA=1.25) and imaged by a CCD camera. Depending on 
the lateral position of the convex lens L3 either the back focal plane or the sample plane is imaged. 
BB beam block, LI, 2, 3 convex lenses, f focal length of LI and L2; S, sample plane and image 



thereof; F, F' back focal plane and image thereof 



image not only the direct space information but also the momentum corresponding to the 
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4l| that such a 



Fourier space. It is indeed a well known fact of Fourier optics [l?} H . 
mapping of the wave vector distribution (as given by Eq. 20) is in principle always possible 
by recording the LR light in the back focal plane F of the oil immersion objective. In 



the optical setup shown in Fig. 9 [39j we realized a dual microscope able to image SPP 
propagation in both the direct and Fourier space. In particular the back focal plane F of 
the oil immersion objective O2 imaged onto a CCD camera in Fig. 9. With such a microscope 




FIG. 10: (A) Direct space LRM imaging of two SPPs beams launched from the ridge and prop- 
agating to the lower left (L) and upper right to be partly reflected (R) and transmitted (T) by 
the Bragg mirror. The image is taken with no beam block. Ridge and Bragg mirror are indicated 
by the dashed line and rectangle, respectively. (B) Fourier plane LRM imaging corresponding to 
(A). C indicates the wave vector components of the directly transmitted laser beam C; R,L,T 
have the same meaning as in (A). (C)-(D) are Fourier filtered images after removing the central 
beam C with a beam block in F 1 (see text). In (C) the reflected beam R is removed by blocking 
the arc of the SPP ring R shown in (B). In (D) the transmitted beam T and the beam incident 
on the Bragg mirror are removed by blocking the arc of SPP ring T in (B). Data from [39] (©, 
Appl. Phys. Lett., American Institute of physics, 2006). 



it is furthermore possible to act experimentally in the Fourier space image plane F' . First 
we can thereby remove the directly transmitted laser beam by using a beam block located 
on the optical axis. As an application of this method of filtering we consider the reflection of 
a SPP beam by an in-plane Bragg mirror. SPPs are launched as before from a gold ridge (50 
nm height, 150 nm width) lithographed on a 50 nm thick gold film. The Bragg mirror 41 1 
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constitutes a one-dimensional lattice of parallel gold ridges (50 nm height, 140 nm width) 
separated by a distance P defining the period of the lattice. The period P is connected to 
the SPP wavelength by Xspp = 2-7T /k S pp < A and to the angle of incidence reflection 9spp 
of the SPP beam relatively to the (in plane) normal to the lattice by 

P = XsP n P . (24) 
2 cos S pp 

In the present case shown in Fig. 10 the Bragg mirror is optimized for A = 800nm (i. e., 
Xspp = 785nm) and for 8spp = 45° incidence angle which means P ~ 555nm. The 
experimental analysis of such a Bragg mirror when the resonance condition (A, Ospp) 
is fulfilled reveals a very high reflectivity of R ~ 95% (see, for example, 34j for some 
earlier results on SPP Bragg mirrors studied with fluorescence microscopy). However, 
in the present experiment we choose an incident angle 9 = 65°. As a consequence the 
reflectivity was much lower (see Fig. 10 A and more details in [39J]). This configuration 




k (1/Mm) 



x (pm) 



FIG. 11: (A) Fourier space cross-cut along the short solid line in Fig. 10 B corresponding to the 
L beam (see inset). Experimental data (symbols) are compared with a Lorentz fit (solid curve). 
The intensity is normalized by setting the maximum intensity of L to 1. (B) Direct space cross-cut 
along the L beam in Fig. 10 A. Data points (symbols) are compared to an exponential fit (solid 
line). Data from [39] (©, Appl. phys. Lett., American Institute of physics, 2006). 

reveals SPP interferences in the vicinity of the mirror (Fig. 10 A). In Fig. 10 B we show the 
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corresponding Fourier space image. The different observed arcs of LR rings correspond to 
the reflected (R), and transmitted plus incident (T) beams. The L beam is associated with 
a SPP launched in the direction to the left, i. e, away from the mirror. C is the directly 
transmitted laser beam distribution. By acting in the Fourier plane image F' of the LRM 
microscope we now block the information associated with the central beam and with the R 
or T beams [39;]. Thereby the according SPP beam images are erased from the image plane 
and consequently weak intensity beams otherwise observed by interference can be directly 
analyzed. For further analysis we extracted radial cross-cuts along the white lines as shown 



in Fig. 10 B [39]. Results are shown in Fig. 11 A for the cross cut along L. The data points 
agree very well with a Lorentz fit given by Eqs. 20, 21. The FWHM of the Lorentzian 
distribution gives us a SPP propagation length of L$pp = 20 /im. This value is in perfect 
agreement with the cross-cut made along the beam L in the direct space image 10 A (see 
Fig. 11 B). Both data agree also with values given by the dispersion relations discussed in 
sect. 1 (see Figs. 2B and 3B). 

LRM was subsequently applied by us to many SPP in-plane devices such as beam splitters 
flol|. dielectric lenses, prisms 18 1 and wave-guides 42 1. In particular, LRM experiments were 
compared to near field optical experiments ( photon scanning tunnelling microscopy) and 



showed good agreement in the cases considered [18|, |42J. LRM appears thus in this context 
as a complementary far-field optical method to NFO such as NSOM. LRM was applied as 
well for analyzing SPP Bragg mirrors (with high reflectivity R ~ 90 — 95%), interferometers 



401 . 143 1 and SPP elliptical cavities 35 1 or 2D SPP microscopes 44 1. In this context we 
observed [40j stationary SPP waves with very high visibility V = (I m ax — -^min) / (^max + Imm) 
by using LRM. This proves directly that SPP wave coherence is conserved in LRM and can 



exploited for quantitative analysis down to the spatial reso 
SPP interferometers such as the ones described in 



34 



35 



ution limit Xspp/2. Actually 



43j reveal clear interference 



pattern and oscillation characteristics of these devices. It is thus possible to develop 2D 
interferometry for SPP waves having all the advantages of current macroscopic interferom- 
etry techniques. 

We also mention the realization of plasmonic crystals (i. e., photonic crystals for SPPs) 
which were studied using LRM (see Fig. 12). In such devices [45| rectangular 2D latices 
made of gold nano-protrusions (200 nm diameter, 50 nm height) deposited on a 50 nm 
thick gold film (see Fig. 12B) are used to create photonic band gaps at specific SPP 
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laser Plasrnonic crystal 



(B 



Au ^ ^ j£/ 




550 nm 





FIG. 12: (A) Principle of SPP in-plane wavelength demultiplexing using plasrnonic crystal. SPPs 
launched from a gold ridge (50 nm height, 150 nm width) are propagating in the direction of a 
2D rectangular lattice (plasrnonic crystal) as described in the text. (B) SEM images of a part of 
the crystal, Bragg lines are indicated by dashed lines. Such crystals can be seen as two sets of 
Bragg mirrors perpendiculary oriented to each other. (C) Sketch representing the unit rectangular 
cell and the direction of the reflected SPP beam for the two different wavelengths at which the 
Bragg reflections occur. (D, E) LRM images of SPP interacting with the demultiplexer at SPP 
wavelengths X x ~ 730 nm and X y ~ 784 nm, respectively. The gold ridge and plasrnonic crystal 
are indicated by the gray line and the gray box , respectively. The vector basis x, y of the latices 



is indicated in A and E. Data from 



45j (© Nanoletters, American Chemical Society, 2007). 



wavelengths X x ~ 730 nm and X y ~ 784 nm (i. e., laser wavelengths of respectively 750 nm 
and 800 nm) corresponding to the two periods of the lattice P x = X x /\/2 = 516 nm and 
P y = X y /V2 = 554 nm. The existence of these band gaps implies that SPP plane waves 
impinging on small devices build up with such lattice will generate stationary waves in the 
crystal. More precisely this implies that SPPs will be reflected in specific and different 
directions if their wavelengths match the values X x or X y and if the angle of incidence 
relatively to the normal to the Bragg planes of the crystal (Figs. 12 B, C) equals 45°. 
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Such devices act consequently as an efficient in-plane wavelength demultiplexer for SPPs 



45] as it was indeed observed experimentally (see Fig. 12D and E). While the results 



discussed here were achieved within the visible spectral range, plasmonic crystal devices 
are expected to perform even better (e.g., in terms of spectral selectivity) in the infrared 
(telecom) spectral range due to significantly lower ohmic losses 46[]. In general, the use 
of multiplexers, splitters and fritters 45) in photonic applications might be specifically 
appealing due to their small footprint in the range of 10 x 10 fim 2 . Furthermore, the use as 



building blocks for classical 



or quantum 



47) optical computing can be envisaged. 



IV. CONCLUSION 

In this article we reviewed the field of leakage radiation microscopy (LRM) theoretically 
and experimentally. Theoretically we analyzed how SPP can generate leaky wave in the 
glass substrate by tunnelling from the air side through a thin metal film supporting SPP 
waves. We showed in particular that for thick film (D > 70 nm) leakage radiation (LR) 
does not affect the dispersion relation on the air/metal interface. Importantly the angular 
distribution of LR is located on a cone matching the SPP dispersion relation. We also 
reviewed the first experimental results reporting the observation of LR on rough surface and 
using near field optics methods. We analyzed more recent application of LRM to SPP nano- 
devices fabricated by electron beam lithography. From all these results we can conclude that 
LRM is a convenient and versatile far field optical method for analyzing quantitatively SPP 
propagation on flat film and their interaction with various nano-devices of direct practical 
interest. Such versatility positions LRM as an appealing alternative to near field optics for 
studying SPP propagation on a scale of, or larger than the wavelength. 

For financial support the European Union, under project FP6 2002-IST- 1-507879 is ac- 
knowledged. 



[1] H. Raether, Surface Plasmons (Springer, Berlin, 1988). 

[2] W. L. Barnes , A. Dereux, T. W. Ebbesen, Nature 424, (2003) 824. 

[3] C. Genet, T. W. Ebbesen, Nature 445, (2007) 39. 



21 



[4] A. Drezet, A. Hohenau, J. R. Krenn, M. Brun, S. Huant, Micron 38, (2007) 427. 
[5] D. W. Pohl, W. Denk, M. Lanz, Appl. Phys. Let. 44, (1994) 651. 

[6] D. Courjon, Near-field microscopy and near-field optics ( Imperial College Press, London, 
2003). 

[7] J. R. Krenn, A. Dereux, J. C. Weeber, E. Bourillot, Y. Lacroute, J. P. Goudonnet, G. Schider, 
W. Gotschy, A. Leitner, F. R. Aussenegg, C. Girard, Phys. Rev. Lett. 82, (1999) 2590. 

[8] B. Hecht, H. Bielefeldt, L. Novotny, Y. Inouye, and D. W. Pohl, Phys. Rev. Lett. 77, (1996) 
1889. 

[9] A. Bouhelier, Th. Huser, H. Tamaru, H. -J. Giintherodt, D. W. Pohl, Fadi I. Baida and D. Van 

Labeke, Phys. Rev. B. 63, (2001) 155404. 
[10] A. Stepanov, J. R. Krenn, H. Ditlbacher, A. Hohenau, A. Drezet, B. Steinberger, A. Leitner, 

and F. Aussenegg, Opt. Lett. 30,(2005) 1524. 
[11] J. J. Burke, G. I. Stegeman, T. Tamir, Phys. Rev. B 33, (1986) 5186. 
[12] J. D. Jackson, Classical Electrodynamics (J. Wiley and Sons, New York, 1975). 
[13] G. N. ZhiZhin, M. A. Moskalova, E. V. Shomina and V. A. Yakovlev in Surface Polaritons, 

edited by V. Agranovich and D. L. Mills (North-Holland Publishing Company, Amsterdam, 

1982) pp. 93-143. 

[14] A. V. Kats, S. Savel'ev, V. A. Yampol'skii, F. Nori, Phys. Rev. Lett. 98, (2007) 073901. 
[15] P. B. Johnson and R. W. Christy, Phys. Rev. B 6, (1972) 4370. 

[16] P. Halevi in Electromagnetic Surface modes edited by A. E. Boardman (John Wiley, New 
York, 1982) Chap. 7. 

[17] L. Novotny, B. Hecht, Principles of Nano- Optics (Cambridge Press, London 2006). 

[18] A. Hohenau, J. R. Krenn, A. L. Stepanov, A. Drezet, H. Ditlbacher, B. Steinberger, A. Leitner, 

F. R. Ausseneg, Opt. Lett. 30, (2005) 893. 
[19] J. A. Dionne, L. A. Sweatlock, H. A. Atwater, A. Polman, Phys. Rev. B 72, (2005) 075405. 
[20] H. J Simon, J. K. Guha, Opt. Comm. 18, (1976) 391. 

[21] C. Sonnichsen, A. C. Duch, G. Steiniger, M. Koch, G. vonPlessen, J. Feldmann, 

Appl. Phys. Lett. 76, (2000) 140. 
[22] M. Brun, A. Drezet, H. Mariette, N. Chevalier, J. C. Woehl, S. Huant, Europhys. Lett. 64, 

(2003) 634. 

[23] M. Brun, S. Huant, J. C. Woehl, J.-F. Motte, L. Marsal, H. Mariette, Solid State Comm. 121, 

22 



(2002) 407. 

[24] H. A. Bethe, Phys. Rev. 66, (1944) 163. 

[25] C. J. Bouwkamp, Philips Res. Rep. 5, (1950) 321. 

[26] C. J. Bouwkamp, Philips Res. Rep. 5, (1950) 401. 

[27] C. Obermmlller and K. Karrai, Appl. Phys. Lett. 67, (1995) 3408. 

[28] C. Obermiiller, K. Karrai, G. Kolb, and G. Abstreiter, Ultramicroscopy 61, (1995) 171. 
[29] A. Drezet, J. C. Woehl, and S. Huant, Europhys. Lett. 54, (2001) 736. 
[30] A. Drezet, S. Huant, and J. C. Woehl, Europhys. Lett. 66, (2004) 41. 
[31] A. Drezet, J. C. Woehl, and S. Huant, Phys. Rev. E65, (2002) 046611. 
[32] T. S0ndergaard, S. I. Bozhevolnyi, Phys. Rev. B69, (2004) 045422. 

[33] H. Ditlbacher, J. R. Krenn, N. Felidj, B. Lamprecht, G. Schider, M. Salerno, A. Leitner, and 

F. R. Aussenegg, Appl. Phys. Lett. 80, (2002) 404. 
[34] H. Ditlbacher, J. R. Krenn, G. Schider, A. Leitner, and F. R. Aussenegg, Appl. Phys. Lett. 81, 

(2002) 1762. 

[35] A. Drezet, A. L. Stepanov, H. Ditlbacher, A. Hohenau, B. Steinberger, F. R. Aussenegg, 

A. Leitner, and J. R. Krenn, Appl. Phys. Lett. 86, (2005) 074104. 
[36] J.-Y. Laluet, E. Devaux, C. Genet and T. W. Ebbesen, J.-C. Weeber and A. Dereux, Opt. 

Express 15, (2007) 3488. 
[37] M. A. McCord and M. J. Rooks, in Handbook of Microlithography, Micromachining and Mi- 
cro fabrication, Volume 1 edited by P. Rai-Choudhury (SPIE and The Institution of Electrical 

Engineers, Bellingham, Washington, 1997), pp. 139-249. 
[38] Saleh, B. E. A., Teich, M. C. Fundamentals of Photonics (Wiley, New York, 1991). 
[39] A. Drezet, A. Hohenau, A. Stepanov, H. Ditlbacher, B. Steinberger, N. Galler, 

F. R. Aussenegg, A. Leitner, J. R. Krenn, Appl. Phys. Lett. 89, (2006) 091117. 
[40] A. Drezet, A. L. Stepanov, A. Hohenau, B. Steinberger, N. Galler, H. Ditlbacher, A. Leitner, 

F. R. Aussenegg, J. R. Krenn, M. U. Gonzalez, J.-C. Weeber, Europhys. Lett. 74, (2006) 693. 
[41] M. Born and E. Wolf, Principles of Optics, seventh (expanded ) edition (Cambridge University 

Press, Cambridge, 1999). 
[42] B. Steinberger, A. Hohenau, H. Ditlbacher, A. L. Stepanov, A. Drezet, F. R. Aussenegg, A. 

Leitner, J. R. Krenn Appl. Phys. Lett. 88, (2006) 09410 . 
[43] A. Drezet, A. Hohenau, A. Stepanov, H. Ditlbacher, B. Steinberger, F. R. Aussenegg, A. 

23 



Leitner, J. R. Krenn , Plasmonics 1, (2006) 141. 
[44] A. Drezet, D. Koller, A. Hohenau, F. R. Aussenegg, A. Leitner, J. R. Krenn , submitted 
(2007). 

[45] A. Drezet, D. Koller, A. Hohenau, J. R. Krenn, Nanoletters (Mai 2007) DOI: 
10.1021/nl070682p. 

[46] T. Nikolajsen, K. Leosson, I. Salakhutdinov, S. I. Bozhevolnyi Appl. Phys. Lett. 82, (2003) 
668. 

[47] E. Knill, R. Laflamme, G. J. Milburn, Nature 409, (2001) 46. 



24 



